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CHAPTER I 


INTRODUCTION 

The TCAS II is a midair threat-alert/coll i sion avoidance system. 
Once every second, the TCAS II system makes a complete survey of other 
aircraft within a range of about 15 miles, measuring their altitudes, 
distance separations, closing rates, and approximate bearings [1]. 

The azimuth bearing of the nearby aircraft is determined by the sum 
and difference patterns of the TCAS II mounted on top of the aircraft. 
Angle calibrations are done by a lookup table in the TCAS II processor. 
This lookup table was created by measurements taken with the array 
mounted on a circular curved ground plane. How well this curved ground 
plane can be used to model the fuselage of an aircraft is crucial to the 
system performance. This report is meant to answer this question. 

As far as the report structure is concerned, it is divided into 
three sections. The first section is used to analyze a circular array 
and to study the scattering properties of a curved and flat ground 
plane. The second section introduces the TCAS II system and the 
corresponding computer model. The last section compares and examines 
the monopulse curves of a TCAS II system mounted on a flat ground plane, 
a curved ground plane and a Boeing 737. Conclusions are made at the end 
of the report. 
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CHAPTER II 


TCAS II CIRCULAR ARRAY ANALYSIS 


A. INTRODUCTION 

Due to the symmetry of a circular array, a beam can be scanned 
through 360° with little change of beamwidth and sidelobe level. The 
ability of a circular array in generating a directional pattern which 
remains constant in shape over a broad bandwidth makes the circular 
array practical for a number of applications. Such applications include 
Direction Finding and Electronic Support Measurement System which is a 
more sophisticated tracking system [2] than the TCAS II studied in 
detail here. 

B. PRINCIPLES OF A CIRCULAR ARRAY 

The far-field patterns of a 2N-element circular array and a 

sin(N<{>) 

2N-element linear array are similar; they both have a . ,, /0 v pattern. 

s i n ^ <p/ c) 

This similarity in far-field patterns means that all the pattern 
synthesis techniques that can be applied to a linear array can also be 
applied to a circular array. The next few sections will discuss and 
make use of this concept. Note that this discussion is based on the 
work of Sinsky, et al . [3] and as such follows his approach to the 
problem. This information is presented here for the convenience of the 
reader. 
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Consider a current ring perpendicular to the X-Y plane as 
illustrated in Figure 2.1. This current ring can be seen as a circular 
array made up of an infinite number of infinitesimal elements. The 
pattern of each of these elements is F(0), and the element spatial phase 
is y(a). The field strength at a point (r,9,<j>) in the far-field is 
proportional to the following integral: 

EU,0}=F(e)a / 1(a) e jY(ct) d a (2.1) 

o 



X 


Figure 2.1. A thin current ring perpendicular to the x-y plane. 
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where 


y{a) = kpa*R = ka(cosax+sinay) *R 
2tra 

or y(a) = “ x~ sine cos(a-4>) . (2.2) 

■■ i k R 1 

Because multiplying constants are unimportant, the factors e J and r - 
are not included. Also, since the excitation function 1(a) for this 
continuous circular array is a periodic function with a period of 2u 
radians, 1(a) can be expressed by a Fourier Series such that [4]: 

N 

I(“) = l C exp(jma) ; N ->• «> . (2.3) 

m=-N m 

Substituting Equation (2.3) into (2.1), one finds that 

EU,0) - F(9)a f J C e J(i"«rtas1necosto-+)) da . ^ 

0 m=-N m 

By the following identity [2]: 

m q 

one obtains that 

N . A 

E ( 4> , e ) = F( 9) 2ua l C j m J (kasine) e JPl4> ; n+<° 

m=-N m m 

Note that J m is the Bessel function of 1st kind, order m. 

When the number of elements in this array is finite, the far-field 
pattern can still be approximated by Equation (2.6). An N-element array 
is drawn in Figure 2.2. It is found that the number of elements in the 


(2.4) 


(2.5) 


( 2 . 6 ) 


4 




Figure 2.2. A N-element array 

array is the same as the number of modes that can be excited [3]. For 
example, if there are eight elements in the array, the number of modes 
will also be eight. In other words, if N=8, the mode number m will be 
equal to -3, -2, -1,0, 1,2, 3, 4, with a total of eight modes. As a result, 
if N is even. Equation (2.3), evaluated at a-j , can be written as 

N/2 2ir 

= l C exp(j— im) (2.7) 

m=-N/2+l 

2ni f h 

where as in Equation (2.3). The excitation of the i element 

can then be written as follows: 
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I. 

1 


( 2 . 8 ) 


N/2 

= l A, exp [j(P + 2 ttU/N)] 

£=-N/2+l X 

where A^amplitude of the £ th mode input 

P = phase of the mode input, and 

A exp(jP) = C in Equation (2.7). 
m r m m 

Based on Equation (2.8), it can be seen that the element 
excitations are determined by the mode inputs. In order to determine 
the correct mode inputs required to produce a desired antenna pattern, 
it is first necessary to remove the mode biases [31. A mode bias is 
defined as the complex radiation intensity (relative to the array 
center) resulting from the excitation of the £ th mode input. Consider 
the approximate far-field pattern given by Equation (2.6) such that 

N m • . 

E ( <|> , 0 ) = F(9)2ira l C m j J m (kasin9)e j * ; N is finite . 

If C = (2iraF(9 ) j m J (kasin9 )) 1 for a particular 9 . then 
m v o m o J o 


E(+) = 


N • A 

y e jm< ^ 

m=-N 


• r 2N+1 . 

sin ( 2~ (j>) 

sin 2 


(2.9) 


where 


(C ) 1 is the mode bias, and C is the mode excitation, 
m m 


This is nearly identical to the formulation for linear arrays. The 
radiation pattern of a linear array of 2N+1 isotropic elements is given 
by 
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( 2 . 10 ) 


2N+1 N 

N sinf - 5 u) 

E(U) = l A e jmu = 2 

m=-N m 


. u 
sin - 


where 


u = (2ua/X)sin<j> , 


4> is the angle off broadside and a is the element spacing of the linear 
array. Comparing the above two radiation patterns, it can be seen that 
the mode compensated circular array will produce the same antenna 
pattern in angle space as that of a linear array in sine angle space. 
Consequently, the use of the Fourier technique for exciting circular 
arrays means that all the pattern synthesis techniques for linear arrays 
may be transferred to circular arrays. In particular, the Tschebyscheff 
linear array weighting function can be applied to calculate the mode 
excitations. 

In conclusion, once the mode biases are known, the element 
excitations can be computed by Equation (2.8). Furthermore, for an 
array with element patterns independent of the angle <f>, the mode biases 
are determined by the elevation angle 9 q of the field point. 


C. THE MODE BIAS CALCULATION 

To compute the mode biases, one needs to know the element 
pattern and the far-field pattern of a N-element array. In general, our 
aircraft code can be used to find the element and far-field patterns of 
an N-element array. However, for simplicity, the example in Sinsky's 
paper [3] is followed as an illustration of the mode bias calculation. 
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Assume the it* 1 element pattern of a circular array is hypothesized by a 
cosine on a pedestal about an axis which is rotated by e s from the 
azimuth angle of the i^ element 4^ and tipped up from the X-Y plane by 
a s [3]. The diagram of this is shown in Figure 2.3. Furthermore, as 
illustrated in Figure 2.4, the element pattern has a shape given by 

F i (Q,+) = l+KX i (2.11) 

where X 1 is the cosine of the angle between the boresight of the i tfl 

-f ->■ 

element (B.. ) and the direction (P). Note that 
K . (10 FB/20 -l) / (10 FB/2n + l) 

with FB defined as the front-to-back ratio. The mathematical expression 
for X.. is given by a simple expression, namely; 



= sine cosa $ [cosU-^-f^)] + cose sina s (2.12) 

If each element pattern has the shape such that 
«s = 0- 8 s ■ 0 


then 


X. = sine [cos(4>-4> i )] 

The far-field pattern from an N-element array is given by 


(2.13) 


B 



(2.14) 


N 

E(e,<j>) = l F (e,0)i exp(jY. ) 

i=l 1 1 1 

where 

F-j ( 9. 4>) is the element pattern of the i th element in the radiation 
direction 

Y.j is the spatial phase of the i^ element relative to the center 
of the circular array which is given by 

Y- = (2na/x) sin9 Ceos ( )] , and 

I.. is the complex element excitation on the i l array element and 
is defined by Equation (2.8). 

After the element pattern and the far-field pattern of an N-element 
array have been computed, the mode biases can also be found. As 
indicated before the mode bias is defined as the far-field radiation 
intensity resulting from excitation of the mode input. For 
convenience, when computing the mode bias, the mode input is defined to 
have a unit amplitude and zero phase. Combining Equations (2.11), 
(2.13), and (2.14), the far-field equation for a circular array, when 0 
is 90°, is given by 

N 2ira 

E( 9=90° ,<f>) = l ((l+Kcosf*-^)) exp j(2irU/N) exp[j— cos($-fj )] 

i=l 

(2.15) 

th 

which is also the mode bias for the a mode input. 
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COORDINATE ( 0,0 = 90°, 


P = sine cos<)> x + sine sin<(> y + cose z 

A 

A A A 

B.. = cosa $ cos(4».+B s ) x + cosa s si n ( <{». +8 $ ) y + sina g z 
Figure 2.3. An element pattern from a circular array. 
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A COSINE PATTERN 


yj«o ° 

I 


* cos y \ 


PATTERN FOR F ( y . ) * I 


* 0 ° 


This far-field pattern may not be realized by any antenna since the 
element patterns defined in Equation (2.11) are only hypothetical. 
However, when the angle 9 is 90°, this far-field pattern gives a good 
approximation to the far-field pattern of a TCAS II array as far as the 
magnitude is concerned. A TCAS II array will be discussed in Chapter 
IV. Table 2.11(a) shows the mode biases for the array with the 
hypothetical element patterns, which give an approximation to the true 
mode biases of a TCAS II array. As has been done in Reference [3], 
these approximate mode biases will be applied to a TCAS II array so as 
to illustrate some basic properties of the circular array. 

If the mode excitations are set proportional to the reciprocal of 
these mode biases, then each mode taken one at a time will produce a 
ejro't’ pattern in the direction (9 0 ,<j>). The overall pattern can be 
approximated by 


sin 


E(<f>) = 




sin(cj>/ 2 ) 


for a particular 0 Q ; 


(2.16) 


whereas, in our example, 0 Q is chosen to be 90°. If the phase 
difference between adjacent modes is <f> o radians, the resulting pattern 
can be expressed as 

• r 2N+1 , m 

sin C 2~(4>-4> 0 )) 

E ( 4 >) = - for a particular 9 0 , (2.17) 

,♦■♦ 0 . 

sin (— 7 -) 

which is the original pattern rotated by <j> Q radians. Note that Equation 

(2.17) can be obtained from Equation (2.18) or (2.19) by changing the 

-jm* 

mode excitation C in Equation (2.20) or (2.21) to C e 
m m 
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D. COMPUTATION OF THE TSCHEBYSCHEFF HEIGHTS 


The use of the Fourier technique for exciting a circular array 
allows the use of all the pattern synthesis techniques for linear 
arrays. To show this explicitly, let us consider the following 
approximate far-field pattern of an N-element circular array: 

N 

N/2-1 j— <*» 

E(4>) = l A m e J ^ + A n/2 e when N is even (2.18) 

m=-N/2+l 


or 


(N-l)/2 

E(*) = l 

m=-(N-l)/2 


A e 
m 


jm<j> 


when N is odd, 


(2.19) 


The above two equations come from Equation (2.9) with Am defined as the 
mode weight, or in this case, the Tschebyscheff weight. When the 
element number (N) is finite. Equation (2.14) gives a much more accurate 
representation for the far-field pattern. However, the approximate 
far-field pattern, expressed by Equations (2.18) and (2.19), is used 
here to show the relationship between the far-field pattern and the 
current distribution. The equation for the current distribution is 
given by 


I ( ) 


N/2-1 jma. 

I A C e 1 
m=-N/2+l m m 


N 

j 2 a i 

+ A C a ^ 

rt N/2 U N/2 6 

when N is even 


( 2 . 20 ) 


or 


I(<*i ) 


(N-l)/2 jma. 

I 

m=-(N-l)/2 m m 


when N is odd 


( 2 . 21 ) 


where C is the mode excitation and A is the mode weight, 
m m 
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Both the current distribution and the far-field pattern are a function 
of the variable (A m ). The manipulation of this variable by pattern 
synthesis techniques can give the desired far-field pattern. To realize 
this far-field pattern, the same variable, or the so-called mode weight, 
needs to be muti plied to the mode excitation to get the correct current 
excitations. Futhermore, when N is even, to obtain a sinN<j>/sin<J> 
pattern, mode N/2 is usually not excited so as to get a symmetry in the 
modes excited. 

If a Tschebyscheff pattern is desired, the mode amplitude 
distribution must be tapered. Consider the far-field pattern of a 
circular array given by 

(N-lJ/2 (N-l)/2 

EU> ’ m= . ( N-l)/2 A ™ 6 = A o + 2 J] V COs(2 *+ /! ° (2.22) 


where A^=A m and N is odd. Since the right side of Equation (2.22) can 
be expressed as a polynomial of order N-l in the variable cos((f>/2), it 
can be matched to the Tschebyscheff polynomial T N ^(U) by letting 

U = U q cos ( 4>/2) 

in which U q is selected so that T^ j(U 0 )=r, with r being the desired 
main beam to si delobe voltage ratio. 

Elliot [5] found that the Tschebyscheff coefficients in Equation 
(2.22) can be approximated by a series such that 


(N-l)/2 

A m = a o + 2 E a cos (p2nm/N) 
p=l v 


(2.23) 
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where 


1 , 1TD. 

a p - Tf T n .j (U 0 cos n ) 
a o ' * T N-1 (U o> 

and is the Tschebyscheff polynomial. 

Combining the above two equations, the mode weights can be expressed as 

1 f (N-l)/2 *p 

A m - N (Vj (U 0 ) + 2 I T N -1 (U o C0S n)c0S (p2irm/N) } (2.24) 

pi 

when N is odd. Similarly, as pointed out in Reference [3], the 
following equation also holds: 

1 N/2-1 *p 

A m = N ^ T N-1 (U o ) + 2 1 T N-1^ U 0 C0S n) cos(p(2m+l)ir/N)} (2 .25) 

p=l 

when N is even. 

E. SUH AND DIFFERENCE PATTERNS 

A circular array of an even number of elements has a radiation 

pattern given by Equation (2.18). Because of the property of 

Tschebyscheff weights, the coefficient (A ) is equal to the coefficient 
, N m 

(A ). Furthermore, since mode 2 is not excited. Equation (2.18) can be 
rewritten as follows: 

N/2-1 

E(*) = A q + 2 l A m cosm<j> . (2.26) 

m=l 

which gives a sum pattern. 
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If current is excited such that 


then 

N/2-1 

E (4>) = A 0 + 2j I A m sinm<t> . (2.27a) 

m=l 

To get a difference pattern, A q must be set to zero. This implies that 
N/2-1 

EU) = 2j l A m sinm<j> . (2.27b) 


The monopulse characteristic equation is given as follows [31: 


M = 201 og 


10 


S+jA 

£-jA 


(2.28) 


By the previous two equations; one can find that 


M = 201 og 


10 


N/2-1 , 

A q + 2 i A m (cosm<|>-sinm<t>) 
m=l 
N/2-1 

A q + 2 i A^ (cosm<t>+sinm<j>) 
m=l 


(2.29) 


F. SIMULATED RESULTS 

A computer program has been written to find the sum and difference 
mode excitations. The results are listed in Table (2.1). Substituting 
these sum and difference mode coefficients into Equation (2.29), a 
theoretical monopulse curves curve is obtained. This curve compares 
favorably with the monopulse curve obtained by UTD from finding the 
far-field pattern of a circular array mounted on a flat ground plane. 
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These two monopulse curves are shown in Figure 2.5. When the target 
angle is outside the -25° to 25° range, the two monpulse curves do not 
agree very well; however, each beam position is only useable by the 
system for a small angular range which is less than ±25°. 

Figures 2.6(a) and (b) show the sum and difference beams when the 
side lobe levels are 50 and 30 dB, respectively. Figure 2.7 shows that 
the sum and difference beams can be rotated without changing the pattern 
shape. 
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TABLE 2.1 


TABLE FOR GENERATING SUM AND DIFFERENCE MODES 


(A) THE MODE BIASES FOR 30 dB SIDE LOBE LEVEL; 
i.e. C” 1= E(90°,0°) from Equation (2.15) 


MODE 

VOLTAGE 
(in volts) 

dB 

DEGREE 

-3 

0.7705 

-2.2649 

-106.3775 

-2 

1.4050 

2.9534 

176.7064 

-1 

1.1957 

1.5527 

123.7852 

0 

0.9377 

-0.5584 

137.7721 

1 

1.1957 

1.5527 

123.7848 

2 

1.4050 

2.9534 

176.7055 

3 

0.7705 

-2.2650 

-106.3784 

4 

0.8570 

-1.3400 

- 33.1587 


(8) THE MODE WEIGHTS; 
i.e. Aj,, from Equation (2.25) 


MODE 

dB 

0 

0.0000 

1 

-1.1716 

2 

-4.9089 

3 

-11.5605 
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TABLE 2.1 (Continued) 


(C) THE SUM MODE EXCITATIONS; 
i.e., f^C m in Equation (2.20) or (2.21) 


SUM MOOE 

dB 

DEGREE 

-3 

-9.2956 

106.3775 

-2 

-7.8623 

-176.7064 

-1 

-2.7243 

-123.7852 

0 

0.5584 

-137.7721 

1 

-2.7244 

-123.7848 

2 

-7.8623 

-176.7055 

3 

-9.2955 

106.3784 

4 

-100. 0000 

33.1587 


(D) THE DIFFERENCE MODE EXCITATIONS; 

.e., A m C m with A^ defined in Equation (2.27a) 


DIFFERENCE MODE 

dB 

DEGREE 

-3 

-9.2956 

286.3874 

-2 

-7.8623 

3.2945 

-1 

-2.7243 

56.2152 

0 

-100.0000 

-137.7721 

1 

-2.7244 

-123.7848 

2 

-7.8623 

-176.7055 

3 

-9/2955 

106.3784 

4 

-100.0000 

33.1587 
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CON 1 COL PATTERN ANGLE 

90 

SLL = 30dB 

BEAM POSITION 

0 



(a) A monopulse curve obtained by UTD from finding the 

far-field pattern of a circular array mounted on a flat 
ground plane. 


CONICAL PATTERN ANGLE 

90 
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0 
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RNGLE 

(b) A monopulse curve obtained by Equation (2.29). 

Figure 2.5. Comparison of monopulse characteristic results. 
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(SCALE: EACH DIVISION- 10DB) 



Figure 2.7. Sum and difference beams with a pattern rotation of 45° 
to the right. 
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CHAPTER III 


HONOPOLE RADIATION FROM A CURVED AND FLAT GROUND PLANE 
A. INTRODUCTION 

The principal analytical tool used to analyze the radiation of an 
antenna mounted on a curved or flat ground plane is the Uniform Theory 
of Diffraction (UTD). Basically, UTD is an extension of Keller's 
Geometrical Theory of Diffraction (GTD) [6]. However, unlike GTD which 
fails in the regions adjacent to the incident and reflection shadow 
boundaries, the UTD solution is valid everywhere. In this chapter, the 
radiation pattern of a TCAS II system mounted on a curved and flat 
ground plane will be analyzed using UTD. The patterns will later be 
compared with those obtained from a TCAS II array mounted on a Boeing 
737 aircraft. The corresponding monopulse characteristic curves will 
also be compared. 

As will be shown later, the TCAS II system can be modeled by a 
circular array of eight pairs of monopoles. According to geometrical 
optics, each monopole will divide the surrounding space into a lit and a 
shadow region by a plane tangent to the surface at the source location. 
Distinct mechanisms are used to calculate the field patterns in these 
two regions. References [7], [8] and [9] provide a good background for 
this material . 
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B. SHADOW REGION 


According to the generalized Fermat's principle, a ray emanating 
from a source, which is located on a curved surface, follows a geodesic 
path on the surface and continually sheds energy into the shadow region. 
Such a creeping wave mechanism is illustrated in Figure 3.1. Point Q is 
the point of diffraction where the ray leaves the surface and propagates 

A 

toward the observation point. As illustrated in Figure 3.2, n' is the 

A A 

unit surface normal vector, t' is the unit tangent vector and b' is the 

AAA 

unit binormal vector given by b'=t‘xn' at the source point O'. Similiar 
vector convention is used along the geodesic path. 

It is found that the field dE radiated by an infinitesimal 
electric current moment dP g at point p g can be expressed as [10] 


dE e (Ps) = d P e (Q')*[n , nT 5 (Q')H+n , b Tg(Q')s] e ‘ jkt e" jks 


~jk z ° / 'Pc / di|o 
t 4ir ' V s( Pc +s) Y dn(Q) 


-Pg(0‘)- 

J^o r T_ 


1 

F 


(3.1) 


where Tg and Tg are the torsion factors at O' and are given in Table 
3.1. H and S are the uniform Fock functions whose expressions are given 
below: 


H = g(e) 

s =i(Q r yg(0 

with 


1 

gU) = / 

®exp(-j2ir/3) 


exp(-jte) 
dx W2'( t ) 


(3.2) 

(3.3) 


(3.4) 


and 
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TABLE 3.1 


PARAMETERS FOR SHADOW RERION 


TYPE or convex 
SURFACE 

SLOT OS if m CAM 

MO NO ROLE OR 
SF # CASE 

SURFACE RAT 
TORSION 

OF CURVATURE IR 
1 DtRCCTlOR 

SURFACE 

diffracted 

RAT CAUSTIC 
OtSTARCt 

K 


T, lO'l 

Tj 10’ 1 

T ( tO'l 

T. (01 

T,(0') 

T.IQ'I 

TIO*» 

^,lCf) 

m 

R 

■ 

0 

B 

B 

0 

0 

• 

• ••»(■$) 

CIRCULAR 

■ 

■ 

■ 

■ 

■ 

jHI 

•lata* 

Q j 


CYLINDER 

■ 

■ 

WES 

■ 

■ 


S« 



CONVEX 

■ 

■ 

T(Q‘l/p <0'l 

H 

1 




t Ji • 

SURFACE 

■ 

■ 

■ 


**/*■ 



Figure 3.1. Perspective view of a surface diffracted ray tube. 
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Figure 3.2. Ray path in the shadow region. 


9 <*> ■ # / 


exp(-jre) 


«°exp(-j2ir/3) 


(3.5) 


which are known as the acoustic hard and soft Fock functions, 
respectively. The Fock type Airy function is given by 


®exp( j2tt/3) 


dt »exp{ rt-t 3 /3) 


I 

and w,, (t) is the derivative of w^ (t) with respect to t. The Fock 
parameter for the shadow region is given by 
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e 


(3.7) 



m(t') 

"^TTT 


with 


m(t‘ ) 


1cpg(t , r| 1/3 

2 


(3.8) 


where pg(t') is the surface radius of curvature along the ray path at 
t'. The width of the surface ray tube at Q is given by 

dn(Q) = P c dip (3.9) 

where p c is the surface diffracted ray caustic distance; t is the length 
of the geodesic path. 


C. LIT REGION 

In the lit region, the radiated energy propagates directly from the 
source to the observation point along the incident ray of geometrical 
optics as shown in Figure 3.3, The electric field at Pj_ in the lit 
region can be expressed as [ 10 ] 

• k - e " JkS 

dE e (P L ) =TJdP e (0')‘ T e *— 5 — (3.10) 

e -jks 

where — § — represents the spherical wave factor associated with the 
point source at O' and dP g has been defined in section R. The explicit 
form of the dyadic transfer function is found to be [ 10 ] 

T* = Z sine 1 [n'n* (HV Vcose 1 ) + n'b T F] . (3.11) 

e o L z o o 1 
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Figure 3.3. Ray path in the lit region. 


The Fock functions H and S for the lit region are defined as follows: 

-je £ /3 -je £ x 

u £ / £. j* e 

H (£ » = -V d£ W ' (T) (3-12) 

and 

-js £ /3 -je £ x 

a z -j e e 

s ( t ) - HTST -77— /_j'2tt/ 3 dT H^M ( 3 . 13 ) 


where the functions and have been defined previously in Equation 

p 

(3.6). The Fock parameter e for the lit region is given by 


e = ” m o (O') cose 1 = -m(O') [1+T cos 0 1 ] 


2 1-1/3 i 

9 COS 9 
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The angle 0 1 is defined by s^n'scose 1 as shown in Figure 3.3. Finally, 
F is given by 


i i i 

S -H cose 

1+T cos e 1 
0 


(3.15) 


as defined in Table 3.2. The other parameters are the same as in the 
shadow case. 


TABLE 3.2 

PARAMETERS FOR LIT REGION 


SLOT OK *; a e«sc 

MONOPOLC OK <? t CASE 

To 

T 

A 

• 

e 

0 

M 

M 

Teoafl' 


B 

T 0 T cm 9 ' 

•U0 [h A ^Tq Teo*® 1 ] 

• la S' T 0 T 

TCO'I Q‘) 

I'-x'c.ff 1 

1 »' 
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D. PATTERN FACTOR 


The general solution for a short monopole mounted on a convex 
surface has been given in the previous two sections. To get a solution 
for a linear antenna problem, one needs to integrate the above solution 
over the source distribution. This gives an accurate solution but is 
very tedious. A more efficient approach is to modify dP^Q') as shown 
in Reference [11] such that 

(a) in the lit region 

P e • n' [1-cos (KL)] (3.I6) 


(b) in the shadow region 


P 

e 


A cos(KLn' •s)-cos(KL) 
n* 7 T~2 

l-(n'«s) 


(3.17) 


where L=the length of the monopole. L should not exceed a quarter 
wavelength for the solution to be valid. 
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E. SOLUTIONS OF A CURVED GROUND PUNE 

The previous three sections gave generalized results to a monopole 
mounted on a convex surface. As illustrated in Figure 3.4, the curved 
ground plane being considered here is part of a circular cylinder. 

Thus, solutions that can be applied to a circular cylinder can also be 
applied to the curved ground plane. Edge diffraction must be included 
in the curved ground plane solution which will be discussed in the next 
section. Combining equations in the previous three sections, Reference 
[121 gives expressions for near-zone fields due to a monopole radiating 
from a circular cylinder as follows 

(a) Lit Region: (see Figure 3.5) 

-jks 

£♦' = P, it ( -(e 1 -t) HV) 1 ♦' (3-18) 


E 0 . - P m { cose 1 • 
-jks 

C Z, Z s e « 

S (e ) s } 0' 


m 


2 1/3“ 

( Q ' ) C k pg ( 0 ' )sin9' ) 


(3.19) 


with 


„ -jZo 

P 1 it = 411 


cos(kln' »s)-cos(kL) 
[l-(n's) 2 ] sin(kL) 


(3.20a) 


and 


o kasin8\i/3 * * 

z = -( 2 ) (n'-t) . 


(3.20b) 
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Figure 3.4. The curved ground plane is part of the surface of a 
circular cylinder. 



Figure 3.5. Ray geometry in lit region. 
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(b) Shadow Region: (see Figure 3.6) 


E = P . . { T. g(e.) 

n dark 1 i=1 ^ 


PgCO^ 


1/6 -jk(£ i +S i ) 
e ' 


/S i (S.+£ i ) n i t (3.21a) 


2 


E b ■ p dark ( & 


-j c 


2 1 / 3 ' 


k Pg (0' )cosW i 


g ( e j ) 


Pg(Oi) 

pJToT 


i/6 -jk(t.+S.) 

e ' 1 


b. 


✓S i (S i +Jt i ) i 


(3.21b) 


with 


P dark 


-jZo l-cos (kL) 

4 it i sinfkL) i 


(3.22) 


Note that the expressions in the lit region are derived by 
employing Ivanov's result [12], Thus, these expressions are not the 
same as those obtained by UTD. However both expressions agree with the 
measured results as shown by Wang [12). 
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Figure 3.6. Ray geometry in shadow region. 

F. EDGE DIFFRACTION FROM A CURVED GROUND PLANE 

The Uniform Geometric Theory of Diffraction developed by 
Kouyoumjian and Pathak [71 was used to solve the edge diffraction 
problem. UTD is a ray optical technique which allows one to gain 
insight into the physical mechanisms involved in the scattering of the 
structure. As shown in Figure 3.7, the incident ray will diffract as a 
cone of rays satisfying the law of diffraction; namely, one finds that 

A A A A 

I# e = d*e (3.23) 
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Figure 3.7. A cone of diffracted rays from a curved ground plane. 


where I is the incident unit vector from the source to the diffraction 
point, 

A 

d is the diffraction unit vector, and 

A 

e is the unit vector parallel to the diffracting edge. 

An edge fixed coordinate system is introduced to define the 
polarizations of the incident and diffracted fields. This coordinate 
system can be expressed as follows: 

A A A A A 

<j>' = -e x I / |e x I| 

A I A A 

8 0 = ♦' X I 

AAA A A 

<J> = e x d / |exd|, and 

A A A 

3 0 = ♦ x d (3.24) 
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As illustrated in Figure 3.2, the surface diffracted field dE e at 
point p s in the shadow region can be expressed as follows [14]: 

-jkZ 0 


dE (P ) = . 

e x s' 4 


KZ 0 - ■* * * A .il/r / PC 

— dP e (Q')-[n'nH 0 ♦ n'bT 0 S 0 le 


where 


^ ^ h / dtlO h 

H o = I L p (Q') exp[-jkt - / c£(t')dt'] 'YdnTQy °p(Q) and 


and 


s 0 = j L P <(! ' ) "Pt-i 1 * - (° «p(t’)< lt 'lV^f5T D p (0) 


djo 


(3.25) 


s,h s,h 

Note that Dp and ap are the Keller's surface diffraction and 

attenuation coefficients respectively [6], They are both related to the 

s ,h 

Airy function, and L p are the lauching coefficients. Equations (3.1) 
and (3.25) are two asymptotic representations to the surface diffracted 
field in the shadow region. In the deep shadow region. Equation (3.25) 
works better; in the region near the shadow boundary. Equation (3.1) 
gives a better result. The surface ray field between Q' and Q is not 
the true field on the surface; it is a boundary layer field [61. If the 
surface ray field comes to an edge, the true field of this incident ray 
can be approximated by Equation (3.25) with the modification of the 
following parameters [14]: 

D p(Q) = 1 D® dge Lp(Q) 


and 


(3.26) 


edge | 

s 1 3D S | s 1 

Dp(Q) = 2 [ ]|(j)'=0 ^p^^ sinfJ Q * (3.27) 
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s n 

Again, Dp and Dp are the diffraction coefficients for the soft and hard 

edge edge 

cases, respectively; whereas, D s and 0^ are the usual edge 

r b s 

diffraction coefficients [15]. Furthermore, L p and L p are the 
attachment coefficients. By reciprocity, they are the same as the 
launching coefficients. 

The expression for the incident ray can be simplified if the 
following relationship is utilized [141: 


1 

2 


I ^(Q')e 
P ^ 



.ctpCndf 


1 1 

p 


Lp(Q' )e 


-/q‘ ap(t ' )dt ' 


Lp(Q) = v(e)//t, and 


(3.28) 


(3.29) 


where v(e) and u(e) are the generalized Fock integrals defined in 
Reference [15]. 

The resulting diffracted field can be expressed as follows: 

-d kZo - 

dEe(Ps) = - j TT dPe(Q')- 


-n''(n‘'4')^ d9e (*,f=0)r(e) 

-n‘ , (b.^)5 o T 0 4 rD s d9e (+.+' =0 ) 


-ju(e) 

ktsinSo 


,-jkt 

7f~ 


[<U> o~* [ 

V 


p 

(p+s)s 


■jks 


(3.30) 


Note that p is the second caustic location related to the edge 
diffracted field and is defined in Equation (3.47). All other 
parameters are the same as those in the shadow region unless they have 
been specified otherwise in this section. 
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Expressions for the soft and hard edge diffraction coefficients are 
given as follows: 

_ e -jir/4 _ 

D s,h = 2n/2iHT sinl^ COt ^ 2n ) F ( KL £ a ^ 

+ cot (TFT) F(KL £ a“(D) * { cot C 2 ^) F(KL £ a + (g + )) • 

+ cot (“^r) F(KL £ a (0 + )) } (3.31) 

where the minus sign applies for D g and the plus sign applies for the 
D^. The angle parameter 0 is given by 

0 + = <j> + ♦' 
and 

+ + 2irnN ± -0 ± 

a'(e') = 2 cos* 2 (3.32) 

The number (N) is desined as the integer which most nearly satisfy the 
following equations [16]: 

2irnN + -0 _ = tt 
and 

2irnN - 0~ = -ir 

Then, 0 _ is associated with the incident field and 0 + is associated with 
the reflected field. The wedge angle number (N) is given by 
WA 

n = 2 " tt (3.33) 
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and the wedge transition function is expressed as [16] 

• 00 rs 

F(x) = 2j /x e JX / e‘ jt dt (3.34) 

/x 

The magnitude and phase of this transition function is shown in Figure 
3.8. The distance parameter (L ) which is equal to tsin 2 a for the 

& V 

curved ground plane case. Also other parameters are defined as follows: 

A A 

-I d*n 
<t> = ir - tan' 1 — 
d»t 

<J>‘ = 0 

_i * * 

= cos (e*d) 
and 

n = 2. 

6. SOLUTIONS OF A FLAT GROUND PLANE 

Figure 3.9 shows the structure of a flat ground plane. The 
solution for a monopole mounted on the surface is given by: 

(a) Incident Field 

e' jkR 

0 P,.. sin0 — z — 0 < 0 < 90°, and 

M t K 

n 3.35 

0 otherwise 

where is defined in Equation (3.20). 
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MAGNITUDE 



PHASE ( OEGREES) 



(b) Diffracted field 




where 

D^ is defined in Equation (3.31) and 

A 

<t> is defined in Equation (3.24). 


(3.36) 


The second caustic location (p) related to the edge diffracted field and 
is given by 

A A A 

J. 1 n e , (s'-s) 

p " T 1 " " a sin 2 0 Q (3.37) 

A l 

where n is the normal to the edge curvature. The incident field (E . , ) 
e <p 

is given by 

e“ jks ' 

" P 11t “s' — (3.38) 

Finally, the total field is the sum of the diffracted field and the 
incident field. Some patterns will be run in Chapter V. 
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CHAPTER IV 


TCAS II SYSTEM 


A. INTRODUCTION 

The TCAS II system is a midair collision avoidance system which 
consists of top and bottom mounted, eight -element, electronically 
scanned arrays. These arrays radiate sum and difference beams at 1030 
MHz and are designed to receive replies at 1090 MHz. Transmissions are 
scanned electronically around 360 degrees in azimuth in 5.625 increments 
with a total of 64 different beam positions. The receiver can provide 
target bearing estimates with an accuracy of one degree. The bearing 
estimates are combined with altitude and slant range information to 
generate accurate 3-D target position and velocity estimates [17], 

Figure 4.1 illustrates a TCAS II system in operation [1]. 


B. THE TCAS II SYSTEM 

The TCAS II system employs a pair of 10.5 inch diameter circular 
arrays of eight top loaded monopoles. A Butler Matrix is used to 
generate the sum and difference patterns as shown in Figure 4.2. 
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Figure 4.1. A TCAS system displays a commercial plane (center) and two 
encroaching planes. The transponder on one of the latter 
planes has an altitude-reporting feature (-12 indicates it 
is 1200 feet below the TCAS plane); the second plane's 
transponder lacks this feature, as indicated by the 
question marks [11. 
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Figure 4.2. Ream forming network which uses a Butler matrix to 
determine the element excitations. 

The Butler Matrix operates in such a manner that given the mode inputs, 
the normalized element excitations are given by [31 
1 N/2 

E i = /w I (A.) exp [jP + j2irU/N] (4.1) 

7 i=-N/2+l 1 1 

where 

= complex excitation of the i^ element (volts) 

A = amplitude of the I*"* 1 mode input (volts) 

X 

P = phase of the l tt7 mode input (radians) 

M = complex excitation of the 1^ mode (volts), and 

X 

N = number of elements (in this case N=8). 
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A circular array is shown in Figure 2.2 in which the elements are 
numbered counterclockwise. The x-axis in Figure 2.2 points toward the 
nose of a aircraft when the y-axis points toward the right wing. To 
steer the beam 0 radians in the counterclockwise direction, a negative 
phase gradient is applied across the mode inputs such that 
j[P 0 -ae] 



When the beam is pointing toward the nose of the aircraft, the mode and 
element excitations are given in Table 4.1. The corresponding sum and 
difference patterns, if neglecting the effects of the tail and wings of 
the aircraft, are shown in Figure 4.3. The target reply will be fed 
into a Squinted Beam Monopulse Processor, the output of which gives a 
monopulse curve. The monopulse characteristic curve corresponding to 
the above sum and difference patterns is shown in Figure 4.4. The 64 
monopulse curves corresponding to the 64 beam positions are not quite 
the same curve. To determine the actual azimuth angle of a target, the 
processor output is compared with a lookup table to compensate for the 
differences between the 64 monopulse curves. 

C. COMPUTER MODEL OF THE TCAS II SYSTEM 

The OSU aircraft code is used to simulate the radiation patterns of 
the TCAS II mounted on the top of a Boeing 737. The circular array is 
modeled by eight pairs of monopoles as shown in Figure 4.5. Each 
monopole pair is used to model a top loaded monopole. In order to 
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TABLE 4.1 


MODE INPUTS AND ELEMENT EXCITATIONS FOR I AND A BEAMS 
AT ZERO DEGREE AZIMUTH 

(a) ELEMENT EXCITATIONS FOR E AND A BEAMS 
AT ZERO DEGREES AZIMUTH (TOWARD NOSE OF AIRCRAFT) 


ELEMENT EXCITATIONS: SUM BEAM 


ELEMENT 

MO. 

AMPLITUDE 

PHASE (DEGREES) 

1 


,482 

32.3742 

2 


.109 

37.8493 

3 


.234 

-78.1094 

4 


.232 

-48.9455 

5 


.240 

-78.2374 

6 


.105 

35.5149 

7 


.477 

32.0745 

8 


.592 

9.1814 



ELEMENT EXCITATIONS: DIFFERENCE 

BEAM 

ELEMENT 

NO. 

AMPLITUDE 

PHASE (DEGREES) 

1 


.371 

90.9364 

2 


.537 

135.1673 

3 


.137 

150.721 

4 


.104 ' 

4.03708 

5 


.259 

-.270649 

6 


.543 

-42.3221 

7 


.427 

-88.7265 

8 


.021 

-30.3929 
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TABLE 4.1 (continued) 


(b) MODE INPUTS TO BUTLER MATRIX TO PRODUCE THE Z 
AND A BEAM ELEMENT EXCITATIONS 


MODE INPUTS: SUM BEAM 


MODE NO. 

AMPLITUDE (VOLTS) 

PHASE (DEGREES) 

-3 

.152 

-99.9442 

-2 

H 

CM 

• 

' -19.8993 

-1 

.473 

46.14 

0 

.644 

3 

1 

.478 

45.163 

2 

.214 

-20.545 

3 

.152 

-99.671 

4 

.00101 

0 


MODE INPUTS: DIFFERENCE BEAM 


MODE NO. 

AMPLITUDE (VOLTS) 

PHASE (DEGREES) 

-3 

.209 

76.0642 

-2 

.254 

151.01 

-1 

.683 

-142.86 

0 

.102 

4 

1 

.537 

39.163 

2 

.324 

-24.545 

3 

.161 

-105.67 

4 

.00101 

0 
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Mil 

««u* wen aivitiw-itMi 


NORMALIZED TO 8.0 dB 
ETH POLARIZATION 
FREQ. - 1.06 GHz 
THC - 90.0 
PHC - 0.0 
THETA - 75.0 


(a) Sum beam 


■MI 



nu 

ncaiti tacN h»i*immoo*> 


(b) Difference beam 


Figure 4.3. Sum and difference beams pointed in the nose direction. 
The aircraft is modeled by a fuselage only. 
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(a) Ideal TCAS II monopulse characteristic curve. 



(b) Expanded monopulse characteristic curve showing 
linearity around antenna boresight. 


Figure 4.4. Ideal TCAS II monopulse characteristic curve and its 
expanded version. 
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TOP VIEW 


NOSE 



ARRAY ELEMENT TAIL 

SIDE VIEW 



Figure 4.5. Geometry of the computer model that was used to simulate 
the circular array. 
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include the coupling effects between the top loaded elements, the model 
has been designed as follows. First, the two monopoles in every pair 
are separated by X/4 and out of phase by 90°. Second, the radiation 
pattern from each monopole pair will give a 15 dB front to back ratio. 
The above requirements can be achieved by exciting the inner monopole of 
each pair with a signal of amplitude (.698) and phase (45°) and the 
outer monopole with amplitude (1.0) and phase (-45°), as defined by 
Sinsky [18], 

To generate a pattern cut, the aircraft code is executed eight 
times. Each execution gives a pattern cut due to one pair of elements. 
The sum and difference patterns are obtained by first multipling the 
corresponding element excitation to each element pattern, then summing 
the eight weighted patterns to produce a steered j or 4 beam. 
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CHAPTER V 


ANALYSIS OF TCAS II SYSTEM ON A CURVED SURFACE 
A. INTRODUCTION 

The objective of this chapter is to answer the following questions: 
Can a curved ground plane be used to model the structure of an airplane; 
or, at least, can it be used to model the fuselage of an aircraft? To 
answer these questions, two programs have been written. The first 
program simulates the far-field patterns of a TCAS II mounted on a flat 
ground plane. The second program models the TCAS II mounted on a curved 
ground plane. The simulated structures are shown in Figure 5.1.* The 
simulated results will be compared with those computed by the aircraft 
code. The aircraft used for this study was the Boeing 737 as shown in 
Figure 5.2. The computer model of the Boeing 737 and its input data set 
are shown in Figure 5.3 and Table 5.1, respectively. Furthermore, it is 
also necessary to define the geometry of the simulated patterns as shown 
in Figure 5.4. The 737 was chosen because its structure is relatively 
simple, and the simulated results agree well with the measured data. 


*Note: For the convience of the reader, all figures and tables referred 

to in Chapter V have been grouped together at the end of the 
chapter. 
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B. PATTERNS FROM A CURVED GROUND PLANE 


Before making any comparisons, the success of the curved ground 
plane program must be verified. With an infinite radius of curvature, 
the curved ground plane results should be the same as those of a flat 
ground plane. For this purpose, the element patterns and the monopulse 
characteristic curves for the TCAS II mounted on a curved and flat 
ground are shown in Figures 5.5 and 5.6. The corresponding sum and 
difference patterns are shown in Figures 5.7 and 5.8. Comparing these 
patterns, one should observe that the curved ground plane patterns agree 
very well with those for the flat one. Note that the flat ground plane 
results have been verified by comparisons with previously published 
results. 

If the sum and difference beams of a TCAS II on a flat ground plane 
are rotated clockwise by 90 degrees, the monopulse curve in each case is 
found to be unchanged as shown in Figure 5.9. Thus, it can be concluded 
that for the same elevation angle, the Butler Matrix will not affect the 
monopulse characteristic curves even though the beams are steered in 
different directions. If the monopulse curves change when the beam 
direction changes, it must be due to the curvature of the surface. 

Next, the curved ground plane program is checked with the aircraft 
code when the surface radius of curvature is finite. The line-of-sight 
signal from the TCAS II array mounted on a curved ground plane should be 
similar to the line-of-sight signal when the TCAS II array is mounted on 
a long ellipsoid as illustrated in Figure 5.10. The input data for this 
ellipsoid is shown in Table 5.2. The sum and difference patterns for 
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each case are il lustrated in Figure 5.11 and the corresponding monopulse 
characteristic curves are shown in Figure 5.12. Comparing the patterns 
between the two cases, it can be seen that the simulated results agree 
well with those using the aircraft code. Note that only the aircraft 
fuselage is present in these results. 

There are still more questions that need to be answered. First, if 
the curvature of a curved ground plane is set equal to the curvature of 
an aircraft at the TCAS II location, then will the curved ground plane 
accurately simulate the system mounted on the fuselage? Next, if the 
diffracted field is also included in the curved ground plane, how well 
does the curved ground plane model simulate the fuselage in terms of the 
resulting patterns? To answer the first question, it is necessary to 
compare the monopulse curves for an aircraft mounted system with that of 
a curved ground plane with the same curvature. First, only the source 
term is used to generate the monopulse characteristic curve of the 
curved ground plane as illustrated in Figure 5.12. It can be seen that 
they agree very well. The corresponding sum and difference patterns 
also agree with each other as shown in Figure 5.11. Thus, it can be 
concluded that if only the source term is considered, a curved ground 
plane can be used to model the fuselage of an aircraft with the same 
curvature. 

To answer the second question, the total field of a curved ground 
plane, that is the sum of the incident and edge diffracted fields, needs 
to be compared with the total field of the fuselage of a Boeing 737, 
that is the incident and surface diffracted fields. Figure 5.13 shows 
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the element patterns for the two cases. The sum and difference patterns 
and also their corresponding monopulse characteristic curves at 
different beam positions are shown in Figures 5.14 to 5.23. Comparing 
the curves for the curved ground plane versus the Boeing 737 fuselage, 
it can be concluded that the edge diffracted field included in the 
curved ground plane solution causes a slight error in the monopulse 
curves if the curved ground plane is used to model the fuselage of a 
Boeing 737. However, this error is small compared to the error 
introduced by the other structures associated with a Boeing 737 as will 
be shown next. Before leaving this section, it is important to note 
that this curved ground plane solution is valid only for a certain range 
of elevation angles as will be discussed in section D. 

C. PATTERNS OF A BOEING 737 MODEL 

The OSU aircraft code is used to study the radiation patterns of 
the TCAS II system. Ideally, the shape of the the sum and difference 
beams will not change as the array is scanned around in the azimuth 
plane. However, the fuselage, wings, tail and other structures will 
distort the antenna patterns; consequently, the performance of the 
monopulse receiver will be affected. Thus, it is important to study the 
structural effects of an aircraft on the radiation patterns of a TCAS 
II. References [19] and [20] give detailed results for this issue. 

To examine the structural effects of a Boeing 737, Figures 5.24 to 
5.38 show the sum and difference patterns and also the corresponding 
monopulse characteristic curves for the following cases: 
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(a) a complete Boeing 737, 

(b) a 2-ft curved ground plane with a curvature of 6-ft, and 

(c) fuselage only of a Boeing 737. 

Comparing the monopulse characteristic curves for the three cases, it is 
found that (b) and (c) can be used to model (a) moderately well if the 
beam position is not near the tail section. However, when the beam 
positions approach the tail section, the aircraft cannot be modeled by a 
curved ground plane. The most severely affected region is in the 
direction of the tail. The large distortions in the patterns are caused 
by shadowing of the vertical stabilizer. 

There is one more observation to make. How does a monopulse curve 
change with elevation angle? For a fixed beam position, the monopulse 
curves spread to the same extent for a flat or a curved ground plane as 
shown in Figures 5.39 to 5.41. As shown in Figure 5.42, the monopulse 
curve changes quite dramatically for all beam positions and elevation 
angles between 20° and 40°. The change was not due to the elevation 
angle changing based on the results shown in Figures 5.40 to 5.42. The 
major variation of the monopulse curve is due to the beam position as 
verified by the data shown in Figure 5.43. This verifies the need for a 
curved ground plane and the need for a lookup table associated with each 
beam position. 

D. LIMITATIONS OF THE SIMULATION 

This curved ground plane program only works for a particular range 
of elevation angles. At the frequency of 1.09 GHz, the program works 
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only when the elevation angle is between 18° and 75°. Perhaps a 
combination of UTD and some other techniques need to be used for a 
complete simulation. When the elevation angle is larger than 75°, there 
is a sudden increase in the number of diffraction points which causes 
the diffracted field to be discontinuous; i.e., a caustic. When the 
elevation angle is lower than 18°, the diffracted field is blocked by 
the ground plane which once again causes a discontinous diffracted 
field. Also, it should be noted that the curved ground plane solution 
used here is not an exact solution because higher order terms such as 
double diffraction are not included. To verify these results, measured 
data should be taken and compared with the present findings. 
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TABLE 5.1 


INPUT DATA FOR TCAS II 
AIRCRAFT CODE: ELLIPSOID VERSION 


ONrINCHES 

3 

FQs 1.06 G BZ 
1,1.06,1. 

EG: BOEING 737 

77.. 74.. 830.. 308 .56 
T 

0 . ,0 . ,70 . 

GG s ELEMENT 1 

2.64,66.464 

2 

1.39,-45. 

0 . ,0 . ,0. ,2.78,3 

1., -45. 

1.39,135. 

0 . ,0 . ,0 . ,2.78,3 
.698,45. 

PG: RIGHT WING 
«,T 

1.. 75.. 67 .952 

1.. 536. 93 .316 .14 

1.. 536 .93 .379.86 

1.. 75.. 240 .26 
PG: LEFT WINS 
4,T 

1., -75., 240 .26 

1., -536. 93 ,379.86 

1., -536 .93 ,316 .14 

1., -75., 67 .952 
PG: TAIL 

4,T 

77. . 8. 25.618.55 
284.147,8.25,819.056 

284 .147.0. . 683 .696 

77.. 0.. 4 83 .19 
PG: TAIL 

4,T 

77.. 0.. 4 83 .19 

284 .147.0. . 6 83 .696 
284.147,-8.25,819.056 

77., -8.25,618.55 


BO: 

T 

PP: FCLAR ELOT IN IB 
P . 

1,1 .5,3 

PD: AZIMJ7H PLANE 

90.. 0..90. 

0,360,1 

T, 1000000. 

EX: 

SG: ELEMENT 2 
3.73,70. 

2 

1.39.0. 

0.,0.,0. ,2 .78 ,3 

1.. -45. 

1.39,180. 

0 . ,0 . ,0 . ,2 .78,3 
.698,45. 

EX: 

SG: ELEMENT 3 
2.64,73.536 
2 

1.39,45. 

0 . ,0 .,0., 2.78,3 

l.,-45. 

1.39,225. 

0 . ,0 . ,0 . ,2 .7 8 ,3 
•698,45. 

EX: 

SG: ELEMENT 4 
0.,75. 

2 

1.39,90. 

0 . ,0 . ,0 . ,2 .78,3 

l.,-45. 

1.39,270. 

0 . ,0 . ,0 . ,2 .78 ,3 
.698,45. 


EX: 

SG: ELEMENT 5 
-2.64,73.536 
2 

1.39,135. 

0 . ,0 . ,0 . ,2 .78,3 

l.,-45. 

1.39,-45. 

0 . ,0 . ,0 . ,2.78,3 
.698,45. 

EX: 

SG: ELEMENT 6 
-3.73,70. 

2 

1.39,180. 

0 . ,0. ,0., 2.78,3 

1., -45. 

1.39.0, 

0. ,0 . ,0 . ,2.78,3 
.698,45. 

EX: 

SG: ELEMENT 7 

-2.64,66.464 

2 

1.39,225. 

0 . ,0 . ,0 . ,2 .78,3 

1., -45. 

1.39,45. 

0 . ,0 . ,0 . ,2.78,3 
•698,45. 

EX: 

SG: ELEMENT 8 
0.,65. 

2 

1.39,270. 
0«,0.,0.,2 .78,3 
l.,-45. 

1.39,90. 
0.,0.,0. ,2.7 8 ,3 
.698,45. 

EX: 
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NOSE 


• s A MONOPOLE 


ELEMENT #1 *8 



(a) A TCAS II mounted on a 2-foot radius curved ground 
plane with a curvature of 6-feet. 


NOSE 



(b) A TCAS II mounted on a flat ground plane. 


Figure 5.1. TCAS II mounted on a curved and a flat ground plane. 
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Azimuth Conical Pattern Angle = 9^ 
Elevation Angle = 90°-9p 
Beam Position = a 


Figure 5.4. Coordinate system for azimuth conical pattern cuts. 
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HOSE 



Radius of disc = 24" 
Elevation angle = 20° 
Frequency = 1.09 GHz 


Element #8 


(a) Element patterns. 



(b) Monopulse characteristics. 


Figure 5.5. Comparison of element patterns and monopulse 

characteristics for a curved and a flat ground plane. 
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NOSE 






Hi 





Radius of disc = 24" 
Elevation angle = 20° 
Frequency = 1.09 GHz 


Element #2 


TAIL 


(a) Element patterns. 



Figure 5.6. Comparison of element patterns and monopulse 

characteristic for a curved and a flat ground plane. 
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NOSE 


HOSE 



TOIL TOIL 


Figure 5.7. The sum and difference patterns corresponding to 
Figure 5.5. 



Figure 5.8. The sum and difference patterns corresponding to 
Figure 5.6. 
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Figure 5.9. Same monopulse characteristic curve when beam positions 
are 0°, 90°, 180°, and 270° for a flat ground plane. 


TCAS I LOCATION 



Figure 5.10. The simulated structure of a long ellipsoid. The source 
field of this structure should be close to the source 
field of a curved disc with a curvature of 72”. 
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TABLE 5.2 


INPUT DATA FOR A LONG ELLIPSOID 


UN: INCHES 
3 

FQ: 1.09 GHZ 
1,1.09,1. 

FG: 

72.. 72.. 1000. .1000. 
T 

0 . , 0 . , 0 . 

BO: 

T 

PP: POLAR ILOT IN IB 
F 

1,1 .5 ,3 

PD: AZIMUTH ILANE 
90 . ,0 • ,70. 

0,360,1 

T.1000000. 

S3: ELEMENT 1 

2.813.-3.536 

2 

1.39,-45. 

0. ,0. ,0. ,2.78,3 

1.. -45. 

1.39,135. 

0 . ,0 . ,0 . ,2 .7 8,3 
.698,45. 

EX: 

SG: ELEMENT 2 

3.979.0. 

2 

1.39.0. 

0 . ,0 . ,0 . ,2 .7 8 ,3 

1 . ,-45 . 

1.39,180. 

0 . ,0 . ,0 . ,2 .7 8 ,3 
.698,45. 

EX: 


SG: ELEMENT 3 
2.813,3.536 
2 

1.39,45. 

0 . ,0 . ,0 . .2 .78 ,3 

1. ,-45. 
1.39,225. 

0 . ,0 .,0. ,2.78,3 
.698,45. 

EX: 

SG: ELEMNT 4 
0.,5. 

2 

1.39,90. 
0.,0.,0. ,2.78,3 

1., -45. 
1.39,270. 

0 * .0 • ,0 • ,2 .7 8,3 
.698,45. 

EX: 

SG: ELEMENT 5 

-2.813,3.536 

2 

1.39,135. 

0 • , 0 . , 0 • , 2 . 78 ,3 

1 . , -45 . 
1.39,-45. 
0.,0.,0. ,2.78,3 
.698,45. 

EX: 


SG: ELEMENT 6 
-3.979.0. 

2 

1.39,180. 

0 . ,0 . ,0 . ,2 .7 8,3 

1. , -45. 

1.39.0, 

0 . ,0 . ,0 . ,2 .7 8,3 
.698,45. 

EX: 

SG: ELEMENT 7 
-2.813,-3.536 
2 

1.39,225. 
0.r0.,0. ,2 .7 8 ,3 

1. , -45. 

1.39,45. 

0 . ,0 . ,0 . ,2 .7 8,3 
.698,45. 

EX: 

SG: ELEMENT 8 
0. ,-5. 

2 

1.39,270. 

0 . ,0 . ,0 . ,2 .7 8,3 

1 . ,-45 . 

1.39,90. 

0 . ,0 . ,0 . ,2 .7 8 ,3 
.698,45. 

EX: 
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(a) Sum patterns 


Frequency = 1.09 GHz 


NOSE 



TAIL 


Source field from a 

curved disc with a 
curvature of 72" 

Source field from 

a long ellopsoid 

. Source field from 
the fuselage of a 
Boeing 737 


(b) Difference patterns 


Figure 5.11. Comparison of sum and difference patterns. 
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(b) Monopulse characteristic curves when beam position is 
90°. 


Figure 5.12. Comparison of monopulse characteristic curves referring 
to the three different cases in the last figure. 
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(a) Element #8 

NOSE 



TOIL 



(b) Element #1 


Element pattern of 
a curved disc 

Element pattern of 
the fuselage of a 
Boeing 737 


Frequency = 1.09 GHz 
Elevation = 20° 


(c) Element #2 


Figure 5.13. Element patterns of a curved disc and the fuselage of a 
Boeing 737. 
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(a) Sum patterns 

NOSE 



TAIL 


Curved disc 

Fuselage of 

a Boeing 737 

Frequency =1.09 GHz 
Elevation = 20° 


(b) Difference patterns 


Figure 5.14. Sum and difference patterns of a curved disc and the 
fuselage of a Boeing 737 when the beam position is 0°. 
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Frequency = 1.09 GHz 


(a) Fuselage of a Boeing 737 

(b) 2-foot radius curved disc 
with a curvature of 6-feet 



Figure 5.15. The monopulse curve corresponding to Figure 5.14. 
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(a) Sum patterns 



(b) Difference patterns 


Curved disc 

Fuselage of 

a Boeing 737 

Frequency = 1.09 GHz 
Elevation = 20° 


Figure 5.16. Sum and difference patterns of a curved disc and the 

fuselage of a Boeing 737 when the beam position is 45°. 
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Frequency = 1.09 GHz 


(a) Fuselage of a Boeing 737 

(b) 2-foot radius curved disc 
with a curvature of 6-feet 



Figure 5.17. The monopulse curve corresponding to Figure 5.16. 
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(a) Sum patterns 
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TAIL 


Curved disc 

Fuselage of 

a Boeing 737 

Frequency = 1.09 GHz 
Elevation = 20° 


(b) Difference patterns 


Figure 5.18. Sum and difference patterns of a curved disc and the 

fuselage of a Boeing 737 when the beam position is 90°. 
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- 16 . -12 


Frequency = 1.09 GHz 


(a) Fuselage of a Boeing 737 

(b) 2-foot radius curved disc 
with a curvature of 6-feet 



Figure 5.19. The monopulse curve corresponding to Figure 5.18. 
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(a) Sum patterns 


NOSE 



TAIL 


Curved disc 

Fuselage of 

a Boeing 737 

Frequency = 1.09 GHz 
Elevation = 20° 


(b) Difference patterns 


Figure 5.20. Sum and difference patterns of a curved disc and the 

fuselage of a Boeing 737 when the beam position is 135°. 
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Frequency = 1.09 GHz 


(a) Fuselage of a Boeing 737 

(b) 2-foot radius curved disc 
with a curvature of 6- feet 



Figure 5.21. The monopulse curve corresponding to Figure 5.20. 
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(a) Sum patterns 


Curved disc 

Fuselage of 
a Boeing 737 



Figure 5 . 22 . Sum and difference patterns of a curved disc and the 

fuselage of a Boeing 737 when the beam position is 180°. 
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Frequency = 1.09 GHz 


(a) Fuselage of a Boeing 737 

(b) 2-foot radius curved disc 
with a curvature of 6-feet 



Figure 5.23. The monopulse curve corresponding to Figure 5.22 



Frequency = 1.06 GHz 
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(a) a Boeing 737 (b) a curved ground plane 
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TAIL 


(c) the fuselage of a Boeing 737 


Figure 5.24. Sum patterns of a Boeing 737, a curved ground plane and 
the fuselage of a Boeing 737. 
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Frequency = 1.06 GHz 


(SCALE: EACH D I V I S 1 0N = 1 ODB) 


NOSE NOSE 
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(a) a Boeing 737 


(b) a curved ground plane 


NOSE 



TAIL 


(c) the fuselage of a Boeing 737 


Figure 5.25. Difference patterns of a Boeing 737, a curved ground plane 
and the fuselage of a Boeing 737. 
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Frequency = 1.06 GHz 



ANGLE 


(a) Boeing 737 

(b) 2-foot radius curved disc 
with a curvature of 6-feet 

(c) Fuselage of a Boeing 737 
Figure 5.26. The corresponding monopulse curves. 
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Frequency = 1.06 GHz 
(SCALE: EACH DIVISI0N-10DB) 
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(a) a Boeing 737 00 a curved ground plane 


NOSE 



TAIL 


(c) the fuselage of a Boeing 737 


Figure 5.28. 


fference patterns of a Boeing 737, 
d the fuselage of a Boeing 737. 


a curved ground plane 
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Frequency = 1.06 GHz 



(a) Boeing 737 

(b) 2-foot radius curved disc 
with a curvature of 6-feet 

(c) Fuselage of a Boeing 737 


Figure 5.29. The corresponding monopulse curves. 
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Frequency = 1.06 GHz 
(SCALE: EACH 0IVISIQN-100B) 
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(a) a Boeing 737 


(b) a curved ground plane 
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TAIL 




(c) the fuselage of a Boeing 737 


Figure 5.30. Sum patterns of a Boeing 737, a curved ground plane and 
the fuselage of a Boeing 737. 
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Frequency = 1.06 GHz 


(SCALE: EACH 0 1 V I S 1 0N= 1 OOB) 

NOSE NOSE 
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(a) a Boeing 737 (b) a curved ground plane 
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TAIL 


(c) the fuselage of a Boeing 737 


Figure 5.31. Difference patterns of a Boeing 737, a curved ground plane 
and the fuselage of a Boeing 737. 
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Frequency = 1.06 GHz 



ANGLE 

(a) Boeing 737 

(b) 2-foot radius curved disc 
with a curvature of 6-feet 

(c) Fuselage of a Boeing 737 
Figure 5.32. The corresponding monopulse curves. 
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Frequency = 1.06 GHz 
ISCflLE: EACH DIVISI0N«1OOB) 


NOSE NOSE 
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(a) a Boeing 737 (b) a curved ground plane 


NOSE 



TAIL 

(c) the fuselage of a Boeing 737 


Figure 5.33. Sum patterns of a Boeing 737, a curved ground plane and 
the fuselage of a Boeing 737. 
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Frequency = l * 06 6Hz 


(SCALE: 


EACH 


oivision»ioo 0) 


NOSE 



TAIL 

(a) a Boeing 737 


(b) a curved ground plane 
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( C ) the fuselage of a Boeing 737 


Figure 5.34. 


Difference patterns 
and the fuselage of 


of a Boeing 737, 
a Boeing 737. 


a curved ground plane 
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Frequency = 1.06 GHz 



(a) A Boeing 737 

(b) 2-foot radius curved disc 
with a curvature of 6-feet 

(c) Fuselage of a Boeing 737 


Figure 5.35. The corresponding monopulse curves. 
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Frequency = 1.06 GHz 


ISCALE; EACH 0IVISI0N-10DB) 
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(a) a Boeing 737 (b) a curved ground plane 
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(c) the fuselage of a Boeing 737 


Figure 5.36. Sum patterns of a Boeing 737, a curved ground plane and 
the fuselage of a Boeing 737. 
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Frequency = 1.06 GHz 
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(a) a Boeing 737 


(b) a curved ground plane 
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(c) the fuselage of a Boeing 737 



Figure 5.37. Difference patterns of a Boeing 737, a curved ground plane 
and the fuselage of a Boeing 737. 
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Frequency = 1.06 GHz 



(a) Boeing 737 

(b) 2-foot radius curved disc 
with a curvature of 6-feet 

(c) Fuselage of a Boeing 737 


Figure 5.38. The corresponding monopulse curves. 
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Figure 5.39. All the monopulse curves of a flat ground plane 

when the elevation angle varies from 20° to 40° and 
the beam position changes. 
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Figure 5.40. The monopulse curves of a curved ground plane with a 
curvature of 6-feet when the beam position is 90° and 
the elevation angle varies from 20° to 40°. 
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Figure 5.42. The monopulse curves of a curved ground plane with a 
curvature of 6-feet when the elevation angle varies 
from 20° to 40° and the beam position changes. 
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Figure 5.43. The monopulse curves of a curved ground plane with a 
curvature of 6-feet when the elevation angle is 30° 
and the beam position changes. 
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CHAPTER VI 


CONCLUSIONS 

A curved ground plane can be used to model the fuselage of an 
aircraft with a slight error in the monopulse characteristic curves. 

This error is small compared to the error caused by the appendages found 
on an aircraft. This is particularly true when the beam position is 
near the tail of an aircraft where the monopulse characteristic results 
are severely distorted by the vertical stabilizer. To reduce such an 
undesirable scattering effect. Reference [20] points out that one 
solution is to move the TCAS II location as far from the vertical 
stabilizer as possible. 

Some other observations made in the previous chapter include the 
followings: 

(a) For a flat ground plane, the monopulse characteristic curves 
do not depend on beam positions due to the symmetry of the 
structure. 

(b) For a curved ground plane, the surface curvature affects the 
monopulse characteristic curves such that they vary 
significantly with beam position. This verifies the need for 
a lookup table for each azimuth beam position. 

(c) For a fixed beam position, the monopulse characteristic curves 
do not change significantly with elevation angle provided that 
the azimuth look angle for each beam is kept small; i.e., less 
than ±25°. 
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Based on the above observations, it can be concluded that no matter 
where the TCAS II is placed on the fuselage, variations in the monopulse 
characteristic curves due to elevation angle changes are rather weak in 
terms of the fuselage curvature, but they can be dependent on the 
aircraft structure such as the vertical stabilizer. The choice of a 
flatter surface does not seem to improve the above condition. However, 
since major errors come from the scattering of the tail, the best 
location of the TCAS II is still to move it as far from the vertical 
stabilizer as possible as suggested by Grandchamp [20j. 

As a last remark, the curved ground plane can be considered as a 
good simulation model for the fuselage of an aircraft but may not be 
good enough to model a whole aircraft. Finally, it should be noted that 
all the conclusions made above are only valid when the elevation angles 
are within the range of the curved ground plane solution. Thus, 
measured results are needed to verify these conclusions. 
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